Baltimore, Maryland 

NOISE-CON 2010 

2010 April 19-21 


Impact of damping uncertainty on SEA model response variance 

Noah Schiller^ 

Randolph Cabell b ) 

Structural Acoustics Branch 
NASA Langley Research Center 
Hampton, VA 23681 USA 

Ferdinand Grosveld c ) 

Lockheed Martin 

NASA Langley Research Center 

Hampton, VA 23681 USA 

Statistical Energy Analysis (SEA) is commonly used to predict high-frequency 
vibroacoustic levels. This statistical approach provides the mean response over an ensemble 
of random subsystems that share the same gross system properties such as density, size, 
and damping. Recently, techniques have been developed to predict the ensemble variance 
as well as the mean response. However these techniques do not account for uncertainties in 
the system properties. In the present paper uncertainty in the damping loss factor is 
propagated through SEA to obtain more realistic prediction bounds that account for both 
ensemble and damping variance. The analysis is performed on a floor-equipped cylindrical 
test article that resembles an aircraft fuselage. Realistic bounds on the damping loss factor 
are determined from measurements acquired on the sidewall of the test article. The 
analysis demonstrates that uncertainties in damping have the potential to significantly 
impact the mean and variance of the predicted response. 

1 INTRODUCTION 

Statistical Energy Analysis (SEA) is often used to predict the high-frequency vibroacoustic 
response of aerospace vehicles. This energy based method requires significantly fewer degrees of 
freedom than conventional displacement based approaches (e.g. FEA); however predictions are 
only available in terms of mean response levels. Mean response levels can be computed over 
various populations, including averages over a population of similar structures, spatial averages 
over a subsystem, and frequency averages within a frequency band. While it is useful to know 
the mean response, in some cases it is important to know bounds on the response. For instance, 
the mean response is rarely used for design since actual response levels can deviate significantly 
from the mean. Therefore when design loads or test criteria are needed, additional factors are 
commonly added to the mean SEA value to account for ensemble, point-to-point, frequency, 
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and/or structural parameter variations. These factors are often based on experience or empirical 
databases 1 . 

Instead of adding margins to the output, there has been significant interest in calculating 
response bounds within the framework of SEA. Much of this work has focused on calculating 
the ensemble variance of built-up systems 2,3 . Recall that SEA predicts the mean response over a 
population of similar structures, such as the mean response within similar vehicles leaving an 
assembly line. Therefore measurements acquired on a particular test structure will not 
necessarily match predictions. Ensemble variance quantifies the expected variability within a 
population of similar structures. 

Lyon and DeJong presented one of the first comprehensive formulations for the ensemble 
variance of built-up systems. While instructive, subsequent work has shown that an underlying 
assumption regarding the statistical distribution of resonance frequencies is not appropriate for 
many systems 4,5 . Langley and Cotoni 3 have since updated and expanded the ensemble variance 
derivation. Their analysis begins with the traditional SEA equations, however the coefficients in 
the equation are taken to be random and a perturbation expansion is used to determine the 
variance of the subsystem energies. VA One 6 software has incorporated /DQJOH\DDQGD &RWRQL^VD 
formulation for ensemble variance. Therefore confidence intervals, accounting for ensemble 
variation, can now be predicted with commercially available software. 

While ensemble variance is valuable, it accounts for just one source of uncertainty. 

Variations in input parameters, such as damping, are not included within the ensemble variance 
formulation. Although Langley and Cotoni outlined a procedure for propagating damping 
uncertainty through an SEA model, the importance of such an approach was not explicitly 
demonstrated 2,7 . The present paper shows that damping loss factor uncertainty can significantly 
impact both the mean and variance of the predicted response. The analysis is performed on a 
stiffened, floor-equipped, composite cylinder that resembles an aircraft fuselage. The paper 
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begins with a brief description of the test structure, which was recently tested at NASA Langley . 

The measurement procedure used to estimate the damping loss factor from experimental data is 
then described. The next section describes the procedure for combining ensemble values, as 
calculated by VA One, with loss factor uncertainty. Finally an example is provided that 
demonstrates the importance of propagating loss factor uncertainty through SEA. 

2 LOSS FACTOR UNCERTAINTY 
2.1 Test Structure 

Loss factor estimates were acquired on the sidewall of a floor-equipped composite cylinder. 

The test structure, shown in Fig. 1, is 1.676 m in diameter and 3.658 m long. The 0.00170 m 
thick shell, which is made of carbon fiber filaments embedded in an epoxy resin, is reinforced by 
orthogonal supports (i.e. ring frames and stringers) that are riveted and bonded to the skin. These 
stiffeners segment the cylinder into 198 rectangular panels with nominal dimensions of 0.203 m 
by 0.381 m. A 0.0127 m thick plywood floor is installed 0.544 m above the bottom of the 
cylinder and the entire structure is supported by two 0.0889 m thick particle board end caps. 

More details on this structure can be found in Jackson 9 . 





Fig. 1 - Photograph (left) and sketch (right) of the stiffened composite cylinder. 


2.2 Impulse response decay method (IRDM) 


Total loss factor estimates for the sidewall were acquired using the impulse response decay 
method (IRDM) 10,11 . This method is used to estimate damping from frequency response 
functions (FRFs). In this case, thirty-six unique FRFs were acquired with drive and 
measurement points spread over two different panels, as shown in Fig. 1. Specifically the 
frequency responses were between a force gage, attached to a shaker on the interior of the 
sidewall, and laser vibrometer measurements acquired on the outside of the cylinder. 

Figure 2 illustrates the steps involved in IRDM. Each frequency response function (shown 
with blue in the figure on the left) is windowed to obtain the band-limited response (shown by 
the red curve). The window is shaped to match the magnitude response of a 4 th -order 
Butterworth band-pass filter with cutoff frequencies aligned with the edges of the 1/3 -octave 
band of interest (shown with the dashed lines). The inverse Fourier transform is then used to 
compute the impulse response (shown with the blue line in the center figure) from each 1/3- 
octave band windowed-FRF. The normalized Schroeder decay curve (red curve on the right) is 
then calculated from the impulse response 12 . A linear regression is used to determine the 
parameters of a simple exponential model that best fit the decay curve 

$<) = c,<r 2 ' / ’", (l) 

where a- 1 is the estimate of the decay curve, rj is the loss factor, and / is the center frequency of 
the band . The square of the correlation coefficient, R , is then computed for each regression to 
quantify the correlation between the original curve and the fit. In this case a perfect fit is 
indicated by R =1. Since most curves are more complicated than a simple exponential decay, the 
results of the curve-fit are sensitive to the fitted portion of the decay curve. Therefore a 
constrained optimization is used in each 1/3 -octave band to select the fit range that gives the 
highest R 2 value. During the optimization, the fit range is constrained between -1 and -20 dB of 
decay. Once the optimal fit range is determined and the fit is performed, screening is used to 
eliminate frequency responses with an R 2 value of less 0.98. The purpose of this step is to 
eliminate bad estimates, which occur when the FRF does not contain sufficient dynamic response 
within the 1/3 -octave band 10 . 





Fig. 2 - Figures highlighting the steps involved in IRDM. The figure on the left shows the 
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the center figure. The figure on the right shows the Schroeder decay curve and the 

curve-fit. 

The mean and variance of the total loss factor were then calculated from all IRDM estimates 
within each 1/3 -octave band from 100 to 4000 Hz. However to calculate confidence intervals it is 
necessary to know the shape of the distribution as well as the mean and variance. Therefore 
histograms of the loss factor estimates were used to estimate the shape of the distribution. 

Figure 3 compares two histograms with normalized log-normal probability density functions, 

where a log-QRUPDO □ GL VWULEXWLRQ □ PHDQ V □ WKDW □ WKH □ YDULDEOH|\ThOR JDULWKP [ 

probability density function (PDF) for a log-normal distribution can be expressed as 
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where p and o are the mean and standard deviation of the variable's natural logarithm. 

Although the number of samples in the plots is too small to definitively characterize the 
distribution of the estimates, it is clear that the data is right-skewed and a log-normal distribution 
appears to be appropriate. 




Total loss factor Total loss factor 

Fig. 3 - Histogram of the total loss factor estimates along with normalized log-normal PDFs in 
the 2500 Hz (left) and 3150 Hz (right) 1/3-octave bands. 

The average loss factor for each 1/3-octave band is shown by the solid blue line in Fig. 4, 
along with the 90% confidence intervals indicated by the dashed lines. Notice that while the 
mean loss factor approaches 0.04 at high frequencies, the 90% confidence interval extends from 
0.015 to 0.083. 



Fig. 4 - The blue curves show the measured total loss factor for the sidewall and the red curves 
show the estimated damping loss factor. The solid lines show the mean values and the 
dashed lines show the 90% confidence intervals. 

The total loss factor estimated by the IRDM includes all energy loss mechanisms, including 
damping and transmission to other components. For example a system composed of two 
identical, coupled subsystems with damping rj and coupling r/ n will have a total loss factor of 

n n + JEn_. ( 3 ) 

This simple example demonstrates that the total loss factor can be significantly larger than the 
damping loss factor if coupling is strong. For complex systems the relationship between the total 
loss factor and damping loss factor is not as simple as Eq. 3. Therefore an SEA analysis was used 
to estimate the damping loss factor from the total loss factor. An SEA model of the cylinder was 
created as described in Appendix A. The damping loss factor was varied within the model while 
estimating the total loss factor as t] T - P/coE , where P is the input power and E is the total 
energy in the driven subsystem. When the estimate of the total loss factor matched the measured 
value, then the damping loss factor was assumed to be correct. This approach assumes that the 
model is accurate. Since the accuracy of SEA tends to deteriorate at low frequencies, predictions 
below 2000 Hz were not used. While the mean values were updated using the SEA model, the 
variance was not. Instead the variance of the damping loss factor was assumed to equal the 
variance of the measured total loss factor. While this approach may not be rigorous, it provides a 
practical way to estimate damping values from total loss factor measurements. Damping loss 
factors are shown with red in Fig. 4. 

The statistics of the damping loss factor can be visualized using the probability density 
function. As an example, Fig. 5 shows the PDF for the damping loss factor in the 2500 Hz 1/3- 
octave band. The circle depicts the mean value while the dashed lines indicate the 90% 
confidence interval. Notice that since the distribution is log-normal, the bounds are not 
symmetric about the mean. 



Fig. 5 - Probability density function for the damping loss factor in the 2500 Hz 1/3-octave 

band. 

3 COMBINED RESPONSE DUE TO ENSEMBLE VARIATION AND LOSS FACTOR 
UNCERTAINTY 

Statistical variation in the loss factor can be combined with ensemble variation using 

h(y) = f f(y\x)g(x)dx, (4) 
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where g(x) is the probability density function for the damping loss factor, and f(y\x) is the 
conditional density of the response 13 . The probability density function for the damping loss 
factor g(x) can be determined from measured data as described in the previous section, and the 
conditional density of the response f(y\x) can be determined using software such as VA One. 

Although VA One provides confidence intervals that account for ensemble variation, the 
conditional density function, f(y\x), is not directly available. As stated in the software 
documentation 6 , the statistical distribution of the subsystem energy is assumed to be log-normal. 

Therefore the confidence intervals are defined as 

(5) 

where E is the response variable (i.e. subsystem energy), p L V □ WKH □ PHDQ □ RI C WKH n YDULDEOH^jV 

logarithm, cr is the standard deviation of the variable's natural logarithm, and q is the 

(1- al 2) -quantile of the standard normal distribution. For 90% confidence intervals, a is 0.1 

and q =1.64. Therefore once the confidence intervals are determined, p and o can be found 

using Eq. 5 and then the conditional density function f(y\x) can be determined with Eq. 2. 

If there are several uncertain system properties, the multivariable extension to Eq. 1 can be 

used, 

h(y) = f f -T /Cy I X l ,x 2 ,...,xjgfx l )g 2 (x 2 )...g n (xjdx l dx 2 ..xbc n , (6) 
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where f(y\xi,X 2 ,...,x„) is the joint conditional distribution of Y given that Xj=xi, X?=X 2 , and X n =x n ; 
and gfxj is the PDF for the uncertain input parameter. By writing the expression in this form, it 
is assumed that Xj, X 2 , through X n are all mutually statistically independent . Unfortunately this 
approach becomes computationally expensive as the number of uncertain input parameters 
increases. This is because it is necessary to resolve the SEA equations and recalculate 
f(y\xi,X20«ryi for every unique combination of (x],X 2 ,...,x„). Therefore in some situations it may 




be helpful to perform an initial sensitivity study to down select to the system properties that have 
the largest effect on the response before proceeding with the analysis. In this paper, the sidewall 
damping loss factor is the only uncertain parameter, so the simple expression given in Eq. 4 can 
be used for the analysis. 

4 RESULTS 

In order to illustrate the impact of damping uncertainty on response level prediction, the 
SEA cylinder model was used to predict panel responses due to a structural power input on the 
sidewall. Specifically the input was applied to panel A, identified in Fig. 9, and responses were 
observed on panels B and C. Panels B and C are located immediately to the right of the driven 
panel as shown in Fig. 9. 

Initially damping uncertainty was neglected and VA One was used to calculate the mean 
response, along with 90% confidence intervals accounting for ensemble variation. In this case 
the mean damping loss factor, shown by the solid red line in Fig. 4, was used for all panels on 
the sidewall. The black curves in Fig. 6 show the mean response (solid line) along with the 90% 
confidence intervals (dashed lines) for panels B and C. The spread in the results is due to 
ensemble variation. 




Fig. 6 - Predicted velocity on panel B (left) and panel C (right). The black curves show the 
output of VA One, which accounts for ensemble variation, the blue curves account for 
damping uncertainty, and the red curves account for both damping uncertainty and 
ensemble variation. The solid lines designate the mean and the dashed lines show the 
90% confidence intervals. (The red curves have been shifted up by 1/3 dB for clarity.) 

Next the impact of damping uncertainty was evaluated without incorporating ensemble 
variation. This was accomplished by solving the SEA equations multiple times while varying the 
sidewall damping loss factor. Again it was assumed that every panel on the sidewall had the 
same damping loss factor and therefore the damping of all panels was changed together. The 
SEA equations were solved one-hundred times with damping loss factors varying from 0.004 to 
0.062. Since ensemble variation was neglected for this analysis, the conditional density functions 
f(y\x ) were represented with Dirac delta functions. The predicted response levels were then 
computed using Eq. 4. The predicted responses for panels B and C are shown with the blue lines 
in Fig. 6. Notice that the spread in the predicted response is much larger when damping 
uncertainty is included in the analysis. It is also interesting to note that the solid blue lines are 
different than the solid black lines. Therefore the mean value of the response changes, by as 


much as 2.5 dB on panel C, when damping uncertainty is included. As a result, it is important to 
account for damping uncertainty within SEA and not just use the mean damping loss factor for 
the analysis. 

Finally the statistical variation in the loss factor was combined with the ensemble variation. 
Once again the SEA equations were solved one-hundred times with damping loss factors varying 
from 0.004 to 0.062. This time the variance package in VA One was used to predict the ensemble 
variation for each run. The 90% confidence intervals predicted by VA One were then used to 
estimate conditional probability distributions f(y\x). Figure 7 shows the conditional density 
functions for three different values of damping x. The three damping values represent the 
minimum, mean, and maximum damping loss factors considered in this analysis. Notice that as 
the sidewall damping increases the peaks shift to the left and become narrower. This indicates 
that both the velocity and the spread in the data are reduced as damping increases. 
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Fig. 7 - Conditional density functions for three different damping loss factors (0.004, 0.018, 
and 0. 062) at 2500Hz in panel B. 

The total response was found by combining the statistics of the damping loss factor, 
characterized by g( x), with the conditional values f(y\x) using Eq. 4. Figure 8 shows the total, or 
unconditional, probability density function of the response, h(y). This PDF accounts for both the 
statistical variation in the damping loss factor and the ensemble variation. The dashed red lines 
show the 90% confidence interval, which was determined by integrating under the curve and 
selecting the bounds such that 5% of the responses were greater than the upper bound and 5% 
were less than the lower bound. 





Fig. 8 - Unconditional density function for the 2500Hz 1 /3-octave band, panel B. 

The combined predictions accounting for both damping and ensemble variation are shown 
with the red lines in Fig. 6. Since the red curves fell almost directly on top of the blue lines, they 
were shifted up by 1/3 dB for clarity. After comparing the blue and red curves, it is clear that 
ensemble variation has little impact on the predicted response. Therefore for this analysis, 
ensemble uncertainty could be neglected without a significant effect on the results. 

4 CONCLUDING REMARKS 

This paper demonstrates the need to account for measurement uncertainty within SEA. The 
analysis was performed on a stiffened, floor-equipped, composite cylinder. The impulse response 
decay method was used to estimate the damping of the sidewall from measured data. Damping 
uncertainty was then propagated through SEA. It was found that uncertainty in the damping loss 
factor has a much larger effect on the predicted response than ensemble variation. The results 
also demonstrate that damping uncertainty can have an impact on both the mean and variance of 
the predicted response. This has important implications. If SEA is performed using only mean 
damping loss factors then the predicted response may be inaccurate. 

5 APPENDIX A: SEA CYLINDER MODEL 

The SEA model of the cylinder is depicted in Fig. 9. The primary components of the model 
are the sidewall, stiffeners, floor, end caps, and acoustic cavities. The sidewall of the cylinder is 
modeled with singly curved shells consisting of 9 layers of 0.0001891 m thick T500-tape in a 
±45, ± 32, 90, +32, +45 layup. The material properties for the tape are shown in Table 1. In 
this case, the 1- and 2-directions lie in the plane of the material, with the 1 -direction in line with 
the fiber orientation. 


Fig. 9 - SEA model of the stiffened composite cylinder without the front end cap or acoustic 

cavities. 


Table 1 - Material properties for T 500-tape. 


Density [kg/m 3 ] 

Elastic Moduli [Pa] 

Poisson Ratios 

Shear Moduli [Pa] 

p 

Ei 

e 2 

e 3 

V'n 

V 23 

^31 

Gl2 

G 23 

G31 

1590 

1.51el 1 

8.96e9 

8.96e9 

0.3 

0.15 

0.15 

5.14e9 

5.14e9 

5.14e9 


As depicted in Fig. 9, the sidewall is segmented into rectangular panels by frames in the 
circumferential direction and stringers in the axial direction. These stiffeners are modeled using 
the beams described in Table 2. Although the composite layup was different for the base, web, 
and crown (top) of each stiffener, the beam section definition in VA One only supports isotropic 
materials. Therefore the equivalent properties of the crown were used for the entire stiffener. The 
equivalent isotropic material properties for each type of stiffener are shown in Table 3. Since the 
modal content of the stiffeners is low throughout the frequency range, they make poor SEA 
subsystems 14 . However it was still important to include the stiffeners in the model to account for 
the impedance mismatch between the panels. The energy storage and dissipation within the 
beams was not modeled though. 


Table 2 - Physical properties of the stiffeners. 


Beam 

Moments 

Cross Section 

Shear Center Offset 


Ixx [m 4 ] 

Iyy [m 4 ] 

Jzz [m 4 ] 

Q 22 [m 4 ] 

A\m z ] 

P\m] 

Dx [m] 

Dy[m] 

Stringer 

4.54e-9 

2.20e-8 

2.65e-8 

2.05e-10 

1.37e-4 

0.14 

0 

8.10e-3 

Frame 

7.60e-8 

1.15e-8 

8.75e-8 

2.01e-10 

1.82e-4 

0.21 

-8.77e-3 

-6.95e-3 


Table 3 - Isotropic material properties. 


Material 

Density [kg/m 3 ] 

Tensile Modulus 
[Pal 

Shear Modulus 
[Pal 

Poisson Ratio 

Composite stringer 

1590 

6.03el0 

2.01el0 

0.49 

Composite frame 

1590 

4.88el0 

2.22el0 

0.099 

Plywood 

711 

1.24el0 

4.66e9 

0.33 

MDF 

730 

2.80e9 

1.04e9 

0.35 


The 0.0127 m-thick plywood floor and the 0.0889 m-thick end caps are modeled as SEA 
plate subsystems with a damping loss factor of 0.03. Meanwhile, the acoustic cavities were 
assumed to have a loss factor of 0.01 . 
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